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Abstract. We consider the semiclassical limit of systems of autonomous PDE's in 1+1 space-time dimensions 
in a scattering regime. We assume the matrix valued coefficients are analytic in the space variable and we 
l/^ I further suppose that the corresponding dispersion relation admits real- valued modes only with one-dimensional 

■ polarization subspaces. Hence a BKW-type analysis of the solutions is possible. We typically consider time- 

' dependent solutions to the PDE which are carried asymptotically in the past and as x —> —oo along one 

mode only and determine the piece of the solution that is carried for x ^ +00 along some other mode 
{^JQ. in the future. Because of the assumed non-degeneracy of the modes, such transitions between modes are 

' exponentially small in the semiclassical parameter; this is an expression of the Landau-Zener mechanism. We 

, completely elucidate the space-time properties of the leading term of this exponentially small wave, when the 

semiclassical parameter is small, for large values of x and t, when some avoided crossing of finite width takes 
\l ' place between the involved modes. 

<N ■ 



1. Introduction 



> 

] Various physical models of wave propagation in space and time are modelled by means of linear systems 
of autonomous Partial Differential Equations (PDE's), with smooth or analytic coefficients in the space 
10 [ variable x G ffi". The solutions to such systems are usually difficult to compute in general and one often 

■ resorts to asymptotic studies in the limit where the wavelength involved is short with respect to the typical 
length scale of the problem, on adapted time scales. This regime is often called space-time adiabatic regime 
or semiclassical regime, due to its relevance in Quantum Mechanics. Typical examples of that situation are 

■ the short wavelength approximation of the wave equation, of Maxwell equations and of the Klein-Gordon 
, equation. Similarly, the semiclassical analysis of the Dirac equation, of the Schrodinger equation in solid 

state physics, or for particles with spin in magnetic fields, and the Born-Oppenheimer approximation in 
molecular physics belong to the same type of problems. This is true also for certain Quantum systems 
whose dynamics is constrained in nanotubes or waveguides. Also, the study of shallow water waves in some 
>^ • linearized regime gives rise to the linearized KdV or Boussinesq equations that share similar properties. 
Plasma physics is another source of physically relevant models entering this category. See e.g. [^], |^], |16|, 

[H], 0- 

From the mathematical point of view, it gives rise to singularly perturbed problems for linear systems 
of PDE's. These problems are tackled with success by means pseudo-differential operator techniques and/or 
BKW methods which provide asymptotic solutions up to errors of order 0{e"^), where e is the ratio of length 
scales, and m depends on the peculiarities of the problem. See the monographs [|l^], Q, [10|, |30|, [36| 
for example. 

The first step in the study of autonomous linear systems consists is using separation of variables to 
reduce the problem to a stationary system, parametrized by an energy variable conjugated to the time 
variable. Superpositions of stationary solutions allow to reconstruct solutions to the full time-dependent 
problem. Then, one determines the dispersion relations or modes of the corresponding symbol, and the 
associated polarization subspaces. We will assume that all modes are real valued, i.e. we will consider 
dispersive waves, according to |3^. In the semiclassical limit, when these real valued modes do not exhibit 
crossings as the position and energy parameters vary, the dynamics of the waves decouples inside the 
polarization subspaces in the following sense: to leading order, independent waves driven by the different 
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scalar dispersion relations propagate along the corresponding polarization subspaces without interacting. In 
particular, transitions between isolated modes, or rather between the corresponding polarization subspaces, 
are forbidden in the semiclassical limit. Moreover, in the scattering limit, these semiclassical transitions 
are typically of order 0{e°°), respectively 0{e~^^^) for some F > 0, in a smooth, respectively analytic 
context. See e.g. |Q, fsj], p, ... This phenomenon goes under the name Landau-Zener 

mechanism, according to the analysis of the adiabatic approximation of the time dependent Schrodinger 



equation (in an ODE context) which yields transitions of this order between isolated eigenvalues, ||39[|, PTj , 
[p^], |pO| , 1^,... Let us recall here that in case the modes experience crossings at some point, the transitions 
may be of finite order in e, indeed of zeroth order in some cases, in the semiclassical limit |ll6| . Their 
determination is technically quite different and we do not address these situations. 

Although extremely small, the transitions between isolated modes computed in the scattering limit are 
quite relevant from a physical point of view in the various examples above. It is therefore desirable for an 
ingoing wave prepared at large negative times along one polarization mode, to determine the asymptotics 
as e — > of the part of the waves that propagate for large positive, but finite, times along another mode, be 
it a transmitted or reflected wave. In a semiclassical context, to achieve such a goal one is lead to further 
require the initial wave to be well localized in energy. 

It is the aim of this paper to determine such exponentially small transmitted waves for quite general au- 
tonomous linear systems of PDE's in 1+1 space-time dimensions, when the coefficients are analytic and 
possess limits they reach sufficiently fast as \x\ — > oo. 

While the conditions allowing the determination of exponentially small transitions between isolated 
modes for a variety of physical situations are rather well understood now in a ODE context, or in the 



language and setting sketched above, for stationary solutions, see [22|, [23], |2C], [^], |35|, |14], [^], pl| , 
[PH], [ p^ j pStI , Q, [^, ... , it is well known that the description of inter- mode transitions in a time-dependent 
context requires more work. The only mathematical results we are aware of regarding this issue concern the 
Born-Oppenheimer approximation in molecular physics p!^ ]. 

The paper ||l^ is mainly motivated by molecular physics considerations and the asymptotic descriptions 
provided there rely heavily on peculiarities of the Born-Oppenheimer approximation considered. However, 
as will become clear, the general strategy of the analysis is actually model independent and, at the price of 
sometimes substantial modifications, it can be adapted to fit the various models and situations mentioned 
above. The importance and frequency of the mechanism of inter-mode transitions in various fields of applied 
mathematics is the main motivation for the present work. Our aim is to extract practical conditions on a 
system of PDE's in 1+1 space-time dimensions under which the exponentially small pieces of propagating 
waves describing inter-mode transitions in a scattering regime can actually be computed, in the semiclassical 



limit. In that sense, the present paper can be viewed as a generalization of [19|. 



Let us describe more precisely the autonomous systems we will be dealing with, the type of results we 
get and the underlying strategy we use to prove these results. Since it requires a fair amount of notations 



and hypotheses to give a precise statement of our main result. Theorem |6.l| , we remain at a rather informal 
level in this introduction. 

Let TZ{x,iedt,iedx) be the differential operator: 

(1.1) TZ{x,ieduied^) = ^ Ain{x){ied,)\iedtT , 

Ze{0,...,m},ne{0,...,r} 

where the d x d matrix valued coefficients Ain{x) are independent of t and analytic in x a neighborhood 
of the real axis. Assuming the matrices Ain{x) possess limit as |x| — > oo which they reach fast enough, we 
want to describe the small e behavior of certain solutions 4>{x,t,e) to the evolution equation 

(1.2) Tl{x,iedt,iedx)<t>{x,t,e) = 0, 
for x E R, in the scattering regime t large, in L^(R). 
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The d X d matrix valued symbol, R{x, E, k) corresponding to TZ{x, iedt,i£dx) writes 

(1.3) R{x,E,k)= Yl Ain{x)k^E^, 

«G{0,...,m},nG{0,...,r} 

where we call the dual variables E and k the energy and the momentum variables respectively. The energy 
parameter will be taken in a window A C M specified below. The associated dispersion relations or modes 
are defined as the set of roots {kj{x,E)} of the polynomial equation in k, of degree md, for x G M and 
^ e A, 

(1.4) det R{x,E,k) =0. 

Our main assumption regarding the type of PDE we consider reads as follows: we suppose there exists 
an energy window A such that for all £" G A, and all x S M, there exist md distinct real valued modes 
{kj{x, E)}Q<j<rnd- The associated kernels R{x, E,kj{x, E)), j = 1, • • • ,md are then shown to be one- 
dimensional and their elements, denoted by Lpj{x,E), are the polarization vectors. 

For comparison and illustration purposes, the case considered in [19| corresponds to {iedt + e^d'^jl + 
^oo(a;))0(2;, t, e) = 0, where ^00(2^) = — ^(a^); is the "electronic hamiltonian" , i.e. a, d x d self-adjoint ma- 
trix. It is assumed that V{x) has non-degenerate eigenvalues {61(2;), • • • , ed{x)} with associated eigenvectors 
{(fi{x), ■ ■ ■ , ipd{x)}. For large enough energies E, (|l.4D yields det {E — A;^/2 — V{x)) = 0, which provides the 
real valued modes {— y^2(i? — 61(3;)), • • • , —^J2{E — ed{x)), \Jl{E — ed{x)), • • • , \Jl{E — ed(x))} and corre- 
sponding polarization vectors {}p\{x)^ • • • , (/^^(x), fdi^), • • • , (pi{x)}. 

Our assumption is very close to the definition of linear dispersive systems in nonuniform autonomous 



medium given in |38|, chapter 11. Such linear systems are characterized there by the fact that the dispersion 



relation can be solved in the form of real roots 

E = W{k,x), with dlW{k,x)^0, 



for k real and x € R. This notion is also reminiscent of the strictly hyperbolic equations |37|. In |37|, a 1+1 
first order partial differential equation is called strictly hyperbolic in x if it can be written as: 

d^^ - A{x, t)dt<^ - B{x, t)$ = 0, 

where the matrix A{xq, to) has real and distinct eigenvalues. If A and B only depend on x, these equations 
are of the same type as (|1.2|) for r = m = 1. However, our assumption and this notion are different, in 



general. The author of [^] gives a characterization for strictly hyperbolic systems of the form 



l<m,l+p<m 

in terms of the principal symbol. By contrast, our assumption concerns the total symbol. 
Separation of variables allows to construct solutions to (|1.2| ) by means of the formula 

(1.5) ^{x,t,e):= [ Q{E,e)e-'*^''i:e{x,E) dE, 

J A 

where ips{x, E) is a solution to the energy dependent stationary problem 

(1.6) R{x,E,iedx)^e{x,E) = 0, 
with 

(1.7) R{x,E,i£dx)= Yl AUx)E''{iedJ, 

le{0,...,m},ne{0,...,r} 

and the function Q(-,e) : A ^ C is an energy density which ensures that E belongs to the prescribed 
window A. The dependence of Q{-,e) on the parameter e will be used to localize in energy the waves we 
want to describe. 
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Equation (|l.q ) is a singularly perturbed system of ODE, to which we apply complex BKW techniques. 
Making use of suitably normalized polarization vectors ipj{x, E), we show that the solutions of (1^) can be 
expanded as 



md 



'ipeix,E) = y2cj{x,E,e)e-'^o^^^y'^^'^y/'^j{x,E), 



where the C valued coefficients Cj{-,E,e) satisfy some linear ODE, which we analyze in the semiclassical 
limit e ^ 0. The assumption S E A imply that the factors e~^'fo ^jiy^^)<^y/^ are phases for all x G M with 
distinct kj(x,E), and the coefficients cj are constant in the semi-classical limit, see e.g. |21|, 

(1.9) Cj{x,E,e) = Cj{0,E,e) +0{e), j = l,---,md. 

The hypotheses on the matrices Ain at infinity ensure the existence of the limits (pj{zizoo, E), kj{zizoo, E), 
Cj{zizoo,E,e), and the error term in (|l.9| ) is uniform in x. In particular, the stationary on-shell scattering 
process characterized by the S'-matrix 

/ ci(+oo,£',e) \ 
C2{+oo,E,e) 

S{E,e)c{—oo,E,e)=c{+oo,E,e), where c{+oo,E,e) = 

\Cmd{+oo,E,e)J 

is well defined. Actually, in our analytic framework, the off-diagonal elements of S{E, e) are exponentially 
small, see below. Thus, for |x| large enough, the solutions (|1.5|) of the time dependent equation (1.2) behave 



as 



(1.10) 



{x,t,e) 



md 

j 



I J Q(£;,e)cj(±oo,£;,e)e-^(*-^+^*^^(±~'^))/Vj(±oo,^) dE. 



Assuming the asymptotic dispersion relations E i-^ kj{zizoo, E) are invertible on A, the asymptotic solutions 
( |1.1[1| ) are given by linear combinations of wave packets associated with each mode and corresponding po- 
larization. The property ( |1.9D shows that transitions between modes induced by the evolution are vanishing 
in the semi-classical limit. 



We determine the asymptotics of certain exponentially small transitions between modes for solutions 
that allow to define a scattering process for \x\ large, in a time-dependent set up. Bona fide scattering 
processes require the energy and the modes we are interested to be such that there exists a mode supporting 
ingoing waves on which we start our solution at time — oo and that there exists another mode describing 
outgoing waves at time -|-oo to which transitions are possible. There exist systems of PDE's that intertain 
outgoing solutions or ingoing solutions only. Our results do not provide interesting informations for such 
systems. 



For definiteness, let us assume in that introduction, that the energy of the waves is well localized around 
Eq £ A \ dA and that, for all x G M, dEkj{x, Eq) < 0. Our sign conventions imply that the asymptotic 



group velocities are then positive, see in particular Proposition 7.1. This implies that such waves travel 



from the left to the right and are polarized along (fj{—oo, Eq) in the remote past. Let us further assume the 
mode kn supports outgoing solutions from the left to the right as well, for x ~ -|-cxd. The incoming waves 
are thus characterized for x large and negative by stationary solutions corresponding to Cfc(— oo, E, e) = 5kj- 
Hence, the summand with label n of ( 1.10| ) corresponding to the coefficient c„(-|-oo, E, e), n 7^ j, determines 
the exponentially small piece of the wave for x in a neighborhood of -|-oo that has made the transition from 
mode kj to mode kn in the course of the evolution, for times t large and positive. 



In order to compute the exponentially small asymptotics of the coefficient c„(-|-oo, E^ e), one uses BKW 
techniques. That is one considers the equation satisfied by these coefficients in the complex plane and makes 
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use of their multivaluedness around points of degeneracy of the analytic continuations of certain modes. As 
is well known, the complex BKW method requires the existence of dissipative or canonical domains, e.g. 
[PI, flOtl, |14], lH^, 1^], ... which is not easy to prove. In our setup, we rely on the analysis 



of which proves that in some avoided crossing regime, dissipative domains exist. The notion of avoided 
crossing requires the introduction of another parameter but we don't want to be specific about this regime 
yet. Let us only mention here that dissipative domains exist in particular when the mode kj{-,E) becomes 
almost degenerate with ki{-,E), with either / = j — 1 or / = / + 1, at one point only on the real axis. The 
outcome of the analysis is the asymptotic formula for e ^ 

(1.11) Cn{+oo, E, e) = T(^)e^^(^)/^(1 + 0(e)) 

with S{E) = k{E) + i^{E), j{E) > and t{E) G C*. The exponent S{E) is given by some action integral 
in the complex plane around the relevant complex degeneracy point of the modes kj(-,E) and kn.{-,E), see 



( [4. 91 ), and the prefactor t{E) possesses some geometric meaning [22]. We localize our wave packets in energy 
by considering typically Gaussian energy densities of the form 

(1.12) Q{E,e) = P{E,e)e-^^-^''^'3/(2e)^ 

where P has support in A and diverges at worst like a polynomial in 1/e as e — > 0. 



With these ingredients, we prove in Theorem 6.1 that for t > large enough, and in the L^(M) norm, 
the piece of the wave function that has made the transition from the mode kj to kn, is given in the limit 
e ^ by 

cPnj{x,t,e) ~ e-"*/^iV,(e) / ^-X2ik-k.)V{2e) ^-iitE+ik)+k.+XAk-k*))/e ^j^ 

Jkn(+oo,A) 

(1.13) +0{l/\tf) + o(e3/^e-^'="*/^7V,(e)). 

Here < /3 < 1/2, A; i— > E^{k) is the inverse function of the asymptotic dispersion relation E i— > fc„(+cxD, E). 
The exponent q*, the average momentum k^ and factors Ai, A2 (such that ReA2 > 0) are determined by 
the action integral S and the energy density Q, and the prefactor A^*(e) is polynomial at worse in 1/e. The 
leading term in (1.13) is of positive norm, constant in time, and of order e^/^, up to the prefactors. 



Moreover, as e ^ and |t| — > 00 this wave is essentially carried on a ball centered at x = —dkE^{k^)t of 
radius ^/t, in the sense, see Proposition 7.1. Finally, the error terms are uniform in e and t, respectively. 



The function ( |1.13| ) corresponds to an exponentially small free wave propagating according to the 
dispersion relation E^{k) with Gaussian momentum profile (within the momentum window /c„(+oo,A)) 
centered around /e*. Note that the error terms are negligible only for large enough times, actually expo- 
nentially large times t ~ e^^^, c > 0. Let us emphasize one point revealed by the present analysis and 
that of Ijl^. The average momentum k^ does not coincide with the naive guess ko ~ kn{+(x>,EQ), which 
corresponds to energy conservation. It is actually dependent on the choice of energy density Q. Similarly, 
the exponential decay rate a* is not determined by the function 7 = Im 5 only, but depends explicitly on 
the density Q as well. 

In other words, the piece of the wave function that has made the transition is asymptotically given for 
small e and large times by the solution to the linear evolution equation, in (rescaled) Fourier space, 

(1.14) iedtfit, k) = E+{k)f{t, k), /(O, k) = e""*/^iV,(e)e-^^('=-'=*)'/(2-)e-^i('=-^*)/^Xfc„(+oo,A). 

Finally, we mention that in case E^{k) is quadratic in k, we can further compute the leading term explicitly, 
as in [^], which yields a freely propagating Gaussian, see Lemma Also, our analysis applies to the 
description of exponentially small reflected waves, as will be explained below. 

The rest of the paper is organized as follows. The precise hypotheses on the operator Tl{x, iet,iedx), are 
spelled out in the next Section. Section]^ is devoted to the analysis of the corresponding stationary solutions. 
The BKW method and the avoided crossing situation are presented in Section ^. The construction of time- 
dependent solutions to the original problem and their scattering properties are given in Section |5|. The 
precise semiclassical analysis in the scattering regime of the time dependent asymptotic waves describing 
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inter-mode transitions is provided in Section |6|. Further properties of the asymptotic waves are given in 
Section 0. A technical Section closes the paper. 

2. Hypotheses for the differential operator 

We consider a differential operator defined by (1.1) where a supplementary small parameter 6 is included 
to define the avoided crossing regime in which the gap between certain modes are small: 

(2.1) n{x,iedt,ied,,6)= An(x, <5)(ie9j'(ie90"- 

le{0,...,m},n£{0,...,r} 

We recall that 

V(/,n) G {0,...,m} X {0,...,r}, Vx G M, \/6e[0,do], An(x, 5) G A^d(C) 

and we define: 

R{x,E,k,S)= J2 AUx,5)k'E^. 

«e{0,...,m},ng{0,...,r} 

Now, we describe the hypotheses on the differential operator TZ. 

(HI): There exist y > and do > such that for any 5 G [0, c?o] the matrix valued functions 

z ^ Ain{z, 6) are analytic in a strip pY = {z £ C; |Imz| < Y} and 

{z, 6) 1-^ Ain{z, 6) is for any {z, 6) £ py x [0, do]. 
(H2): There exist v > 1/2, c > and 2mr matrix valued functions {S ^ A;„(iboo, (^)}"^||J^''"^I 

such that V(5 G [0, do],: 

sup \Rez\'^^''\\Ainiz,5) - Aini+oo,6)\\+ sup |Rezp+''P/„(z,5) - A„(-oo,5)|| < c. 

z£pY,Kez<0 zGpy ,Rc2>0 

Now, we describe the assumption of avoided crossing. We assume that A C M is a compact interval 
with non-empty interior such that, for any E € A: 
(H3): For any x G K and any 5 G [0, do]) there are md real values {ki{x, E, 6), k2{x, E,6), . . . , k^^i^, E, 6)} 
such that det R{x, E, k, 5) = 0. 

For any 6 G [0, do], the values kj{x,E,S) have md distinct limits as x — > — oo and as x — > +oo that 
we denote by kj{±oo, E, 6). The labels are chosen as follows. 

When (5 > 0, the functions kj{x,E,5) are distinct for x G [— cxo,+oo] and are labelled by: 

ki{x, E, 6) < /c2(x, E,6) < ■■■ < kmd{x, E, 6). 

When (5 = 0, the functions kj (x, E, 0) such that det R{x, E, k,6) = are given by md real functions 
that have p{E) > finitely many real crossings at xi{E) < • • • < Xp(^E)iE)- Precisely, we assume 
that for some fixed positive Y and for any fixed E £ A that 

• The functions kj{x,E,0) are labelled according to: 

ki{-oo,E,0) < k2{-oo,E,0) < ■ ■ ■ < kmdi- oo , E , 0) . 

• For any j G {1, . . . , md}, the function (z, E) i— > kj{z, E, 0) is continuous on py x A. 

• For any j G {1, . . . , md}, the function z i— > kj{z, E, 0) is analytic on py- 

• For any / G {1, . . . ,p{E)}, there exist exactly two integers G {1, . . . , md}^ such that: 

h{xi{E),E,0) = kj{xi{E),E,0). 

Besides, we assume that 

d,{h-kj){xi{E),E,0) ^0. 
For certain results, we also impose the condition that these avoided crossings be generic in the sense 



of m, 20 



(H4): Fix Eq G A. Near an avoided crossing (xo(-E'o)) -^o) of /cj and kj, there exist three functions 
E a{E), E 1-^ b{E), E c{E) such that, in a neighborhood of Eq, 
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(1) The difference kj — ki satisfies: 

[kjiz, E, 6) - ki{z, E, <5)]2 = a^{E){z - xo{E)f + 2ciE){z - xo{E))6 + b\E)5^ + R^iz - xq{E), 5), 

where i?3 is a remainder of order 3. 

(2) We have: 

a{E) > 0, b{E) > 0, a^{E)b^{E) - c^{E) > 0. 
According to [^], we know a priori that the functions kj are analytic in both variables except at the crossing 
points. The assumptions (HI), (H2) and (H3) imply analyticity in both variables at the real crossing points: 

Lemma 2.1. Assume that (HI), (H2) and (H3) are satisfied. Then, for A small enough, 

(1) the number p{E) does not depend on E ^ IS.. 

(2) there exists y > such that (z, E) kj{z, E, 0) is analytic on py x A, for any j G {1, . . . , md}. 

(3) for I G {1, • • • ,p{E)}, the function E ^ xi{E) is analytic on A. 

We will prove Lemma p.l| in Section |8|. 
Similarly, assumptions (HI) to (H4) imply the following result: 

Lemma 2.2. Under assumptions (HI) to (H4), the functions a, c and 6^ are analytic in a neighborhood of 
Eq. Besides, a{E) = \d^{ki - kj){xo{E), E)\. 

Lemma |2.2| is proven in Section @. 



Let us end this Section by noting here that our hypotheses imply that the modes are real, but they 
do not guarantee that the norm is conserved under the time evolution. This question is addressed in 
Section ^. 

3. Generalized Eigenvectors 

In this Section, we assume that R and A satisfy (H3), and we investigate the properties of the modes, 
their corresponding polarization vectors and the stationary solutions. For the time being, the parameter 
(5 > is fixed and we drop it in the notation. The generalized eigenvectors ips{x,E) G are defined as 
solutions of the time independent equation: 

(3.1) R{x,E,ie^^)^l^e{x,E) = 0. 
For any E G A, the set of such solutions is md-dimensional. 

We define: 

r 

(3.2) V/G{l,...,m} Ni{x,E) = ^Air,ix)E'' 

n=0 

so that 

m 

R{x,E,k) = ^Ni{x,E)kK 

1=0 

We first prove the following result: 

Lemma 3.1. We assume that R and A satisfy (HI) and (H3). We have the following properties. 

(1) For any E G A and any x G M, Nm{x,E) is invertible. 

(2) For j G {1, . . . , m}, (z, E) Nj{z, E) is analytic in py x A. 

(3) // we define H{x, E) by: 

Id ... 

Id ... 



(3.3) H{x,E) 



... Id 

-{N-^No)ix,E) -{N-^Ni)ix,E) -{N-^N^){x, E) ... -{N-^ Nm-i)ix, E) 
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then a{H{x,E)) = {k; det{Y, Ni{x,E)k^) = 0} = {k; det{R{x, E, k)) = 0}. 

1=0 

(4) The functions {{x,E) i-^ -£')}jG{i,.--,md} analytic in M x A. 

(5) Ker {R{x,E,kj{x,E)))) is one- dimensional. 

Proof. The singular values of R{x, E, k) = X^J!^o -^/(^' E)k^ are the roots of the polynomial 

(3.4) L{k) = det A^«(x, . 

This polynomial is of degree md and the highest coefficient is det Nmix, E). According to (H3), since L has 
md distinct roots, det Nm{x,E) ^ 0, which proves (1). 

Assertion (2) is immediate, consider (3). A complex number k G a{H{x, E)) if there exists $ G C™'^\{0} 
such that = Block by block computations show that $ is of the form: 

/ ^ \ 

k{p 

(3.5) $ = . , G C^, 
with det(^™oA^,(x,S)A:') = and 99 G Ker(^[^o iV;(2;, ^)/cO- 

Again by [^], (H3) with 5 > implies that the functions kj{x,E) are analytic in a complex neighborhood 
of M X A which proves (4). Point (5) follows from (|3.5| ) and the fact that a{H{x, E)) is simple 

We introduce some normalized eigenvectors of R{x, E, kj{x, E)). 

3.1. Canonical eigenvectors of R(x, E, kj{x, E)). For a matrix A, we denote its adjoint by A* =* A. 

Fix j G {1, . . . , md}. Under (H3) and according to [^6|, we know that there exist two vector valued functions 
and with values in such that: 

(1) The functions {x, E) ^ E) and (x, E) ^ ^\{^^ ^) analytic on M x A. 

(2) \/{x,E) G Mx A, ij{x,E) G KeiR{x,E,kj{x,E)), i]{x,E) G KeYR*{x,E,kj{x,E)). 
Definition 3.1. The vector (pj = ctjS^j, with 

(3.6) aj{x,E)=e ' <,]M,niu,E,.^H,> 
is called a canonical eigenvector associated to R{x,E,kj{x,E)). 

We notice the following facts: 

• The vector ipj does not depend on the choice of G Ker {R*{x, E, kj{x, E)))). In particular, we can 



choose ft so that 



• Condition ( |3.6| ) may seem artificial but we shall see in the proof of Lemma 3.2 that it corresponds 
to the Kato normalization of the eigenvectors of H{x,E). 

3.2. Decomposition Lemma. 

Lemma 3.2. We assume that R and A satisfy (H3). Let ip^^XjE) be a solution of (p.l|). There exist md 
functions {{z,E,e) 1— > Cj{z, E,e)} j^^i^ ^^j^j .such that: 
(1) The function ijj£{x, E) satisfies: 

md 

V/ G {0, . . . , m - 1}, {ied^yil^six, E) = Y^ Cj{x, E, £)k\{x, E)e--^^o ''i(y'^)'^y^j{x, E). 



(2) If we define 



( Ci{x,E,£) 



c{x, E, e) = 

\ Cmd{x,E,£) 

The vector c satisfies the following differential equation: 
(3.7) dxc{x, E, e) = M{x, E, £)c{x, E, e), 

where the matrix M is given by: 



(3.8) 



(3.9) Vj / I, 



Mji{x,E) = aji{x,E)e' - 

with Aji{x,E)= / [kj{u, E) — ki{u, E)]du, 
Jo 

and \/j G {1, . . . , md}, E) = 0, 

"^'^"^'^^ " kj{x,E)-h{x,E) 
< ip'\,R{x,E,ki)dxipi > +dxki < ip\,[dkR{x,E,ki) - dkRix,E,kj)]ipi > 

< ip],dkR{x,E,kj)ipj > 



(3) For any j G {1, . . . , md}, ipj is a canonical eigenvector of R{x, E, kj{x, Ej) and ip - is any eigenvector 
in KerR*{x,E,kj{x,E)). 

Remark: The set {ipj}j^^i ... md} is a linearly dependent family of vectors in C"'. The decomposition 
in point (1) above corresponds to the familiar BKW Ansatz in semiclassical analysis, see e.g. Q. 
Proof. Let ipi;{x,E) be a solution of (|3.l| ). We define: 



(3.10) 



^six,E) 



{i£dx)Tpe{x,E) 

{ied,)"'-^Mx,E) 



Then ^^{x,E) satisfies: 

(3.11) i£dx^six,E) = H{x,E)-i>e{x,E). 

Equation (|3Tl|) has been studied in [||, H. We use the results obtained there and write: 

md 

(3.12) H[x,E) = ^k,j{x,E)Pj{x,E), 

where the matrix valued functions Pj{x,E) are the one-dimensional eigenprojectors of H{x,E) and satisfy: 

md 

^Pj{x,E)=Imd. 

i=i 

Hypothesis (H3) implies the existence of a basis of eigenvectors of H{x,E): (x, £')}j=i_...^md- 

We determine these eigenvectors uniquely (up to a constant depending on E) by requiring them to satisfy: 

(3.13) H{x,E)<^j{x,E) = kj{x,E)<^j{x,E), Vj = 1, . . . , md, 

(3.14) Pj{x, E)dx<^j{x, E)=0, Vj = 1, . . . , md. 
Indeed, recall that if VF(a;, E) is the solution of 

md 

(3.15) dxW{x,E) = Y,idxPj{x,E))Pj{x,E)W{x,E), W{0,E)=I^d, 



it is well known that W{x,E) satisfies the intertwining identity 

W{x,E)Pj{0,E) = Pj{x,E)W{x,E), Vj G {l,...,md}. 

The generator of ( 3.15| ) being analytic in E, W is analytic in both variables {x,E) G M x A, see section 
XI. 5. Hence, 

^j{x,E) := W{x,E)<^j{0,E), Vj G {l,...,md} 
where {$j(0, E)}j^^i ... is basis of analytic eigenvectors of H{0, E), satisfy 



Pj{x,E)<^j{x,E) = <^j{x,E) and eq. (|3.14D . 



We refer to [|2|, ||, ^ for the details. 
We will rewrite the eigenprojectors as: 

Pj{x,E) 



< <i>Ux,E),<^j{x,E) > 



\^j{x,E) ><<^]{x,E)\, 



where ^]{x,E) G Ker(i?*(x,^) - kj{x,E)), since kj{x,E) = kj{x,E). 

We use the same notation for duality in C"* and C'"'^ since no confusion should arise. 

Let us begin by specifying equation ( p. 14 ) in our case. We consider an eigenvector E) of H{x, E). 
Sj is written as: 



^.j{x,E) 



ij{x,E) 
kj{x,E)Cj{x,E) 

kf-\x,E)C,ix,E) 



, with £,j{x,E) G KerR{x,E,kj{x,E)). 



The vector <I>j must be of the form $j = cxj'Bj, where aj G C and we define tpj = otj^j. Then: 



L J 



Now, if G Ker {H*{x,E) - kj{x,E)) then $j satisfies (^J§ if: 



It remains to choose and to compute < sj, dx'^j > and < sj, >. 

We start with the computation of the vector E). It is an eigenvector of H*{x, E) associated with the 
eigenvalue kj{x,E) = kj{x,E). Let ^Ux,E) G Kei{R*{x, E,kj{x, E))). We check that we can take 



1=1 

m 
1=2 



N* 



Then: 



< -t, - >= ^ k^-' E ^5"' < >= E f^y' E < ^^1' >=< ej, E '^^'^c, > 
p=i /=p «=i p=i 1=1 
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--<^],dkR{x,E,kj{x,E))Cj 



> 



Similarly, we compute: 



< E],d,E, >=< i],dkR{x,E,k,)d.,ij > + < ej, ^dlR{x,E,kj)d^ij > . 

This implies that <pj is a canonical eigenvector of R{x,E,kj). 
From 1 24, we know that any solution to ( |3.11| ) can be written as: 

md 

^e{x,E) =Y,Cj{x,E,e)e-yo''^^y^^'^'^y<^j{x,E), 
where the scalar coefficients cj satisfy the differential equation: dxC = Mc, where M is given by (^.g|), and 



"^■^ " 



We compute 



<<i>],d,'^i> = d.,hY,Y.^q-i)k]-^^"^ <n;^],^i> 

p=2 q=p 
m m 

(3.16) + j2Y.T'^l'<^*M^^^^i>- 

p=l q=p 

By interchanging the indices p and q and according to the formula 

Phi - " " 



a — b 

p+l=s 



we obtain formula (p.S|). The first statement of the Lemma stems from formula ( 3.10| ). This ends the proof 



of Lemma 3.2 



3.3. Behavior of the matrix M. The following Lemma describes the behavior of the coefficients and 
phases entering the definition of M. 

Lemma 3.3. We assume that (HI), {H.2) and (US) are satisfied. Then, 

• The eigenvalues kj satisfy for any A; G N and any I G N: 
(3.17) 

yEeA, sup \x\'^+''\d''Ed^{kj{x,E) -kj{+oo,E))\+ sup \xf +''\d'Ed^{kj{x, E) - kj{- 00, E))\ < 00. 

• The eigenvectors ipj satisfy for any I £ N, uniformly in E £ A: 

(3.18) sup \x\^+''\\d^E{ipj{x,E) - ipj{+oo,E))\\+ sup {xl^+^Wd^Eiipjix, E) - Lpj{-oo, E))\\ < 00. 

• Moreover, for any A; G N* and I £ N, uniformly in E £ A: 

(3.19) sup \x\^+''\\d^Ed^{^j{x,E) - ipj{+oo,E))\\+ sup Ixl^+^Wd^Ed^iy^jix, E) - ^j{-oo, E))\\ < 00. 



• For any A; G N and any I £ N, the coefficients of the matrix M satisfy uniformly in E £ A: 

(3.20) VxGM, V(j,p) G {1,. . . ,md}^ |(9^(9^ajp(x, ^)||xp+'' < 00. 

• Let 

(3.21) ujj{±^,E)= / [kj{y,E)-kj{±oo,E)]dy, 

Jo 
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anc 



(3.22) 



rx 

J kj{y,E)dy = xkj{±oo, E) + ujj{±oo, E) + rf{x, E). 



Then we have, uniformly in E & A, and for any n G N, 



(3.23) Vj G {1,... ,md}, sup\x\'^+''\d%r+ {x, E)\ + sup\x\^+''\d'l;rj {x, E)\ < oo. 



x>Q 



x<0 



We prove Lemma 3.3 in section 8.3. 



3.4. The vector c. In the following lemma, we describe the behavior of the vector c defined by the ODE 

Lemma 3.4. We assume that (HI), (H2), and (H3) are satisfied. 

• For any E £ A and e > 0, the limits Cj{zizoo, E, e) exist for all j = 1, - ■ ■ , md. 

• // the initial conditions to (3.1) are chosen so that c{—oo,E,£) is uniformly bounded in E £ A and 
e > 0, then we have for some constant C uniform in e and E £ A: 

|(9ijcj(±cx),£;,e)| + |cj(±oo,^,e)| < C, 

sup \x\'\dECj{x., E, e) — dECj{+oo, E,e)\ + sup \x\'^\dECj{x, E, e) — dECj{—oo, E,e)\ < C, 

x>0 x<0 

sup |a;|^'^'^|cj (j;, E, e) — Cj{+oo, E,e)\ + sup |rE|"'"^''|cj(x, E, e) — Cj{—oo, E,e)\ < C. 

x>0 x<0 

Remarks: 

i) As the proof shows, the condition sup ||c(— oo, i?, e)|| < oo can be replaced by 



EeA 

3x0 S ^ such that sup ||c(xo, -B, e)|| < oo. 



(3.24) 



ii) In the construction of solutions to (|l.2| ) by means of an energy density, we can (and will) always assume 
that the initial conditions, wherever they are chosen, are uniformly bounded in energy: 



(3.25) 



3xq £ M such that sup ||c(xo, -B, e)|| < oo. 



iii) The equation being linear, we can actually always assume condition ( 3.24| ) holds. This is what we do in 
the rest of the paper. 



We shall prove Lemma 3.4 in section 



According to Lemma 3.4, we can define the stationary scattering matrix S{E,e) by: 
(3.26) S{E, e)c(-oo, E, e) = c(+oo, E, e). 

In order to describe the time-dependent scattering processes we are interested in, we need more detailed 
informations about the stationary ^-matrix. 



4. Complex BKW analysis 

In this section, the parameter (5 > is still kept fixed. All the information about transmissions 
and transitions among the asymptotic eigenstates is contained in the asymptotic values of the coefficients 
Cj{x,E, ±oo) defined in section |3.4| and hence in the stationary scattering matrix S{E,e). We extract 
this information by mimicking the complex BKW method of |24| and |21], while keeping track of the in- 
dependence. 

In the simplest setting, the complex BKW method requires hypotheses on the behavior of the so-called 
Stokes-lines for equation (3.11) in order to provide the required asymptotics. These hypotheses are global in 
nature, and in general are extremely difficult to check. See e.g., However, in the physically relevant 
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situation of avoided crossings, they can be easily checked, as is proven in and will be recalled in the 
next Section. We restrict our attention to these avoided crossing situations. 

To study the 5-matrix, it is enough to consider the coefficients cj that are uniquely defined by the conditions 

(4.1) Cj{—oo,E,e) = l Ck{—oo,E,e)=0, for all k^j. 

The key of the complex BKW method lies in the multivaluedness of the eigenvalues and the eigenvectors of 
the analytic generator H{x, E) in the complex x plane. 

According to (H3), the eigenvalues and eigenvectors of H(x,E) are analytic in x on the real axis. They 
may have branch points in py that are located in 

(4.2) n{E) = {z e PY ; ^ j such that kj{z, E) = ki{z, E)} 

4.1. The set Q{E). By the Schwarz reflection principle, for any £^ G A, we have Q{E) = Q{E). Besides, 
the set y Q{E) is bounded in py- 

EGA 

We have the following description of ^{E), see pl[| : 

Lemma 4.1. Fix Eq £ A. There exists a neighborhood Aq of Eq and a finite number R of bounded open 
sets {f^i}i(={i,...,_R} in py n C+ such that: 

• For any E G Aq, 9.{E) C Uf f^iUf 

• For a/H G {1, . . . , R}, n M = 0. 

• For any E G Aq, and i G {1, • • • contains only one crossing point. This point is a crossing 
point for finitely many distinct couples of modes. 

We define 17 = Uf Uf 

Under our genericity hypotheses, we have the following local behavior at a complex crossing point zq G 
0(Eo),: 

k,{z, Eo) - ki{z, Eo) ~ ^{Eo)iz - zo)'^\l + 0(z - zq)). 

The eigenprojectors of H{x,E) also admit multivalued extensions in py \ ^{E), but they diverge at generic 
eigenvalue crossing points. We only have to deal with generic crossing points. To see what happens to a 
multivalued function f in py\i} when we turn around a crossing point, we adopt the following convention: 
For E fixed, we denote by f{z,E) the analytic continuation of / defined in a neighborhood of the origin 
along some path from to z. Then we perform the analytic continuation of f{z,E) along a negatively 
oriented loop that surrounds only one connected component 0,i of fi. We denote by f{z,E) the function 
we get by coming back to the original point z. We define Co to be a negatively oriented loop, based at the 
origin, that encircles only Qi when fij G C4.. When G C„, we choose Co to be positively oriented. 

We now fix Qi G C+. For any E G Aq, if we analytically continue the set of eigenvalues {kj{z, E)}'^]^^, along 
a negatively oriented loop around ^li, we get the set {kj{z, E)}'^]^^ with 

kj{z,E) = k„g(_j){z,E), for j = !,■■■ ,md, 

where 

(4.3) TTo : {1,2,--- ,m} ^ {1,2,--- ,md} 

is a permutation that depends on Oj. As a consequence, the eigenvectors possess multi- valued analytic 
extensions in py\^. The analytic continuation ^j{z,E) of ^j{z,E) along a negatively oriented loop around 
ilj, must be proportional to <l>7rQ(j)(z, E). Thus, for j = 1, 2, • • • , md, there exists Oj{(o) G C, such that 

(4.4) ^,iz,E) = e-''^^<°'^^<^^,(^,)iz,E). 
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The above imphes a key identity for the analytic extensions of the coefficients Cj{z,E,£), z G Py\^- 
Since the solutions to (|3.1l|) are analytic for all z & py, the coefficients Cj must also be multi- valued. In our 
setting, Lemma 3.1 of implies the following lemma. 

Lemma 4.2. For any j = 1, . . . , md, we have 

(4.5) c,{z,E,e) e-^^co '^.K^)'^-/^ ^"^^^(^0'^) = c^,^,^{z, E, e) 

where Co ^.i^d iTo{j) are defined as above and are independent of E & Aq. 

Remark: Since has a finite number of connected components, it is straightforward to generalize 
the study of the analytic continuations around one crossing point to analytic continuations around several 
crossing points. The loop Co can be rewritten as a concatenation of finitely many individual loops, each 
encircling only one connected component of 0. The permutation ttq is given by the composition of associated 
permutations. The factors e~^^^^'^°'^^ in (|4.4| ) are given by the product of the factors associated with the 

individual loops. The same is true for the factors exp — i /^^ ^ji^, E) dz/e^ in Lemma L2, 



4.2. Dissipative domains. We now describe how to use the above properties in order to control the limit 



e — > 0. The details may be found in [21|. 



The idea is to integrate the integral equation corresponding to ( ^.7| ) along paths that go above (or 
below) one or several crossing points, and then to compare the result with the integration performed along 
the real axis. As z — > — oo in py these paths become parallel to the real axis so that the coefficients take 
the same asymptotic value Cm{—oo,E,e) along the real axis and the integration paths. Since the solutions 
to ( 3.11D are analytic, the results of these integrations must agree as Rez — > oo. Therefore, (4.5) taken at 



z = oo yields the asymptotics of c^Q(^j^{oo, E, e), provided we can control Cj{z,E,e) in the complex plane. 
We argue below that this can be done in the so-called dissipative domains of the complex plane. We do 
not go into the details of these notions because a result of ||2l| will enable us to get sufficient control on 
Cj{z,E,e) in the avoided crossing situation, to which we restrict our attention. 

We recall that Aji is defined in ( p. 81 ). We rewrite ( p. 7] ) as an integral equation: 

(4.6) Cj{x, E,e) = Cj{xo, E, e) + / y^aji{x',E)e ^ ci{x' , E,e)dx' . 

Jxo I 

By explicit computation, we check that ( [4.6| ) can be extended to py \ 0. We integrate by parts in ( |4.6| ) to 
see that ([4.6|) with xq = —oo can be rewritten as: 




(4.7) Cmiz,E,e) = - zey^ .^^^^^^ e^^-(^'^)/^ q(^, ^, 

k[z,E) 

Mf^lJ) ] f^^^^^'^^)/^c,{z',E,e)dz' 

ki{z',E) 

as long as the chosen path of integration does not meet 0. Here, ~ denotes the analytic continuation 
along the chosen path of integration of the corresponding function defined originally on the real axis. This 
distinguishes Crn{oo,E,e) from Cm{(x>,E,e) computed along the real axis as x ^ oo. These quantities may 
differ since the integration path may pass above (or below) points of Q. If the exponential factors in ( [4.7| ) 
are all uniformly bounded when e — > 0, as it is the case when the integration path coincides with the real 
axis, it is straightforward to get bounds of the type 

(4.8) Cmiz,E,e) = + Os(e). 

However, when dealing with Cm in the complex plane, these exponential factors are usually not uniformly 
bounded and one needs to restrict integration paths to certain domains in which useful estimates can be 
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obtained: 



One defines a Dissipative domain for index j, Dj C py \^ associated with the initial condition ([4.1|), 
by the conditions: 

• Dj C py \ ^ and sup^^jj. Rez = oo, infz^Dj Re 2; = —00, 

• For any z G Dj and any index k G {1, . . . , md}, there exists a path 7*^ C parametrized by 
u G (—00, i] which the regularity properties 

lim Re7''(n) = — 00 , ^''{t) = z and sup sup |5„7^(u)| < 00 

^eDj Me(-oo,t] 

• satisfies the monotonicity properties 

u I— > Im Ajyfc(7'^(ii)) is non-decreasing on (— oOjt]. 

Again, as it is well known, the existence of paths from —00 to +00 passing above (or below) points in 
and along which the exponentials can be controlled is difficult to check in general. We can overcome 
these complications by restricting attention to avoided crossing situations where the existence of dissipative 
domains for all indices has been proven in [^], see hypothesis (AC) below. The interest of the definition 
above lies in the following property. 



When a dissipative domains exists for the index j, ( [4.5| ) and (4^) imply 
(4.9) c.„(,-)(oo,£;,e) = e'^^Co ^.K^)'^"/. e-^o,iCo,E) (^i + Oe{s)), 



where the Oe{£) estimate is uniform for E G Aq. This is the main result of Proposition 4.1 in |21] for our 
purpose, under the assumption that a dissipative domain Dj exists. 

In our context, all quantities depend on S G Aq. However, by carefully following the proof of Proposi- 



tion 4.1 of |21], it is not difficult to check that the estimate (|4.8| ) is uniform for E G Aq. For later purposes, 
we also note here that under the same hypotheses on the exponential factors, ^Cm{z,E,e) is uniformly 
bounded for < e < Eq and E G Aq for some fixed eo, by differentiation of ( [4.71 ). See the proof of Lemma 



3.4 for this property on the real axis. 



4.3. Avoided crossings. We now make use of the avoided crossing situation, that allows us to prove the 
existence of dissipative domains. We thus restore the parameter 5 in the notation. We therefore work under 
(H3) and under the following assumption on the patterns of crossings for the modes {kj{x, E,0)}: 

(AC): 

• For all X < xi{E), 

ki{x,E,0) < k2{x,E,0) < ■ ■ ■< k^d{x,E,0). 

• For all J < / G {1, 2, • • • , md}, there exists at most one Xr{E) with 

kj{Xr{E),E,0) - ki{Xr{E),E,0) = 0, 

and if such an Xr{E) exists, we have 

d_ 

dx 



(4.10) —{kj{xr{E),E,0) - ki{xr{E),E,0)) > 0. 



• For all j G {1, 2, • • • , md}, the mode kj{x, E, 0) crosses modes whose indices are all superior to j or 
all inferior to j. 

To any given pattern of real crossings for the group {kj{x, E, 0)}, with E G Aq, we associate a permutation 
TT as follows. The modes {kj{x, E, 0)} are labelled in ascending order at x ~ —00, by (H3). Since there are 
no real crossings for G A and as x — > -|-oo, the values {kj{x, E,0)} are ordered uniformly in G A at 
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X = +00. If kj{+oo, E,0) is the k^^ eigenvalue in ascending order at x = +00, the permutation vr is defined 
by 

(4.11) TT{j) = k. 

Let be in a sufficiently small interval Aq. For a loop Co that surrounds all the complex crossing points 
and ttq the associated permutation (see (^)), ttq corresponds to the permutation vr. 

We can now restate the main result of |2l| that describes the asymptotics of the coefficients defined in ([4.7|). 
We only have to check that, for small 5 > 0, dissipative domains exist and do not depend on £' G Aq. We 
refer to pl| for the details. The construction of these dissipative domains is based on a perturbation of the 



case 6 = 0. By mimicking the arguments of [21|, as in we obtain that estimates of the type ([4.9D are 



true for certain indices j and n, determined by the permutation ( |4.11[ ): 
Theorem 4.1. Assume that (HI) to (H3) are satisfied and that (AC) holds. If 5 > and Aq are small 



enough, the TT{j),j elements of the matrix S{E,e), with 7r(j) defined in (4.11) have small e asymptotics for 
all j = 1, ■ ■ ■ , md given by 



S.U)AE,e) = n e-^'K"^'^)e^4^^'(^'^'^)'^^/^(l + 0,,,,(.)), vr(j) | 



where, for 7r{j) > j (resp. 7r(j) < j), Q, I = j,-- - ,7r{j) - 1 (resp. I = j, • • • ,7r(j) + I), denotes a 
negatively (resp. positively) oriented loop based at the origin which encircles the complex domain il^ (resp. 
^r) corresponding to the avoided crossing between ki{x,E,d) and ki^i{x,E,6) (resp. ki-i{x , E , 6) ) . The 
J^^ ki{z, E,6) dz denotes the integral along Q of the analytic continuation of ki(0,E,6), and 9i{Q,E,5) is the 



corresponding factor defined by m-A) 



Remark: Under our regularity hypotheses in 6, it is easy to get the following properties, see [^ ], 

lim / ki{z,E,6)dz = 0. 

Let us emphasize here that we do not have access to all off-diagonal elements of the S'-matrix; those we 
can asymptotically compute are determined by the pattern of avoided crossings. Moreover, there are cases 
in which one can compute all elements of the S'-matrix, due to supplementary symmetries in the problem, 



see [25|. Sometimes, the coefficients to which we have access are not even the largest ones in the avoided 
crossing situation, as shown in [25|. 

On the basis of steepest descent arguments, transitions between modes that do not display avoided crossings, 
i.e., that are separated by a gap of order 1 as 5 ^ 0, are expected to be exponentially smaller than the 



transitions we control by means of Theorem iA, as 5 shrinks to zero. Since the coefficients in the exponential 
decay rates given by the theorem vanish in the limit 5 — > 0, it is enough to show that the decay rates of the 
exponentially small transitions between well separated levels are independent of 6. 



That is the meaning of the following proposition, which is proven in [19|: 



Proposition 4.1. We assume that (H3) is satisfied. Further assume that the eigenvalues of H{x,E,S) can 
be separated into two distinct groups cri{x, E, 6) and cr2{x, E, 6) that display no avoided crossing for E £ A, 
i.e., such that 

inf dist{ai{x, E,6), a2{x, E, 6)) > g > 0. 

a:£p^U{±oo} 

Let P{x,E,5) and Q{x,E,S) =1 — P{x,E,5) be the projectors onto the spectral subspaces corresponding to 
(7i(x, E, 8) and cr2{x, E, 6) respectively, and let Us{x, xq, E, 6) be the (space-) evolution operator corresponding 
to the equation 

(4.12) ie -^Us(x,xq,E,6) = H(x,E,6) Ue(x,xo,E,5), with Ue(xo,xo,E,6) =1. 
dx 
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Then, for any 5 > 0, there exists £o{5), C{5) > depending on 5, and T > independent of 5, such that for 
all e < eo{6), 

Jim \\P{x,E,6)Ue{x,xo,E,6)Q{xo,E,6)\\ < C((5) e'^/^ 

XQ — > — OO 

This Proposition implies that the stationary transitions between modes without an avoided crossing 
are exponentially smaller than transitions between modes displaying an avoided crossings. It also shows 
that in any case, these transitions are all exponentially small. 

Let us end this section by remarking that we have always specified initial conditions at x = — oo. 
Obviously, the BKW analysis can be equally performed for coefficients whose initial conditions are specified 
at X = +00, mutatis mutandis. 



5. Exact solutions to the time-dependent equation 
In this Section, we construct solutions to 

(5.1) TZ{x, iedt, iedx)4>{x,t, e) = 0, xGM 

by taking time-dependent superpositions of the generalized eigenvectors 'tp£{x,E), for E € A, studied in 
Section |^. We investigate particularly these exact solutions in the scattering regime of large but finite times 
t, and for any fixed e > 0, not necessarily small. 

The superpositions of generalized eigenvectors depend on an energy density Q{E,£) that might be 
complex valued. We assume that the following regularity conditions holds: 

(CO): The density E Q{E,e) is supported on A and is on A, for any fixed e. 
Moreover, (J^lEj ) is true. 

In this Section, the parameter 5 is fixed and we omit it in the notations. We work under the hypotheses 
(HI), (H2), and (H3) and we define: 

„ md 

(5.2) (p{x,t,e)= i;e{x,E)e^Q{E,e)dE = y2(t>j{x,t,e), 
where 

(5.3) 4'j{x,t,e) = / Cj{x,E,e)e ^ ipj{x,E)e ^ Q{E,e)dE. 

J A 



Since the integrand is smooth and A is compact, <f){x,t,e) is an exact solution of ( |1.2| ). 
We also get from the decomposition ( p. 10 ) for alH = 0, • • • , m — 1, 

(5.4) {ied^y(l){x,t,e) = J Cj{x,E,e)e i k'j{x,E)ipj{x,E)e—Q{E,e)dE 



md 

Y.4>f{x,t,e), 
i=i 



with the convention 0^'^' (x, t, e) = 4'j{x,t,e). Note, however, that in general {iedx)^ 4'j{x,t,e) / (x, t, e). 

The behavior of (j)j{x,t,£) for large x can be understood under the following supplementary assumption 
(GV): 

(5.5) Wj £{l,...,md}, V£;gA, dEkj{±oo, E) ^ 0. 

Let us note that condition (GV) is quite natural. Indeed, with our sign conventions, —l/dEkj is the group 
velocity of the asymptotic waves (5.9). Our condition says that we want to describe waves with non-zero 

17 



asymptotic velocity. Moreover, (GV) also imposes the presence of at least one time derivative in the 
definition of the differential operator Tl{x,iedt,iedx, 6). 
We have 

Lemma 5.1. Assume that (HI), (H2), (H3), (CO) and (GV) are satisfied. Let 

^+ = ^^^P lo , /]_ ^ > 

EGA, j£{l,...,md} \OEKj[±00,I^)\ 

and 

^~ = ™f m-TTT > 0- 

EeA, jG{l,...,md} \OEKj{±00, E)\ 

Fix a G (0, 1). Then, there exists > such that, for x large enough and for either t = or any t ^ 
and x satisfying: 

K+ , , , K- 
\x/t\ > or \x/t\ < , 

we have for all j = 1, ■ ■ ■ , md: 

\\4>j{x,t,e)\\<^ and \\^f{x,t,e)\\ < ^, 
\x\ ■' \x\ 

where / € [0, • • • , m — 1] and \\ ■ \\ is the norm in C". 

Specializing to the j'th mode, there exist x^{j) S and Cs{j), independent of time, such that for any 
(5 £ (0, 1) and any I £ [0, • • • ,m — l], if\t\ > 1 with sign (t) = ibsign {dEkj{zizoo, E)) and ibx > ibx^(j) then 

(5.6) Uf{x,t,e)\\ < ^'^^^ 



\tf\x\(^-P)' 
Remarks: 

i) As direct corollaries, we get that (j)^j\-,t,e), and thus {edxY (!){', t,e) belong to L^(M), for any t G M, and 
any / = 0, • • • , m — 1. Moreover, 

(5.7) sup ||(e9,)^,/.(-,t,e)||i2(ffi) =0(C,). 

\t\<i 

ii) The behavior in e of and Cs{j) cannot be estimated under hypothesis (CO) only. However, anticipating 
on our eventual choice of Q{E,e), see ( |6?l| ) below, if the energy density satisfies 

(5.8) sup{\Q{E,e)\ + \edEQ{E,e)\) < oo, 

B6A 
e>0 

and assuming ( p.24| ), then the constants and Cs{j) are actually uniform in e — > 0, as easily checked from 
the proof and Lemma 3^ . 

iii) The complicated looking second statement simply says the following, for < /3 < 1/2. In the asymptotic 
regions where (/>^-' is driven by the asymptotic group velocity —l/dEkj{±oo,E), if time flows in the wrong 
direction, in the sense that the wave is driven out of these regions, then the norm over those regions 
decreases. 



iv) We prove this Lemma in section 8.5 



In a scattering regime, we expect our solutions to behave as freely propagating waves along independent 
modes. Let us introduce such asymptotic waves (p{x,t, e,±oo): 

md 

(j){x, t, e, ±oo) = 4>j{x, t, £, ±oo), 

with 

(5.9) (/>j(x,i,e,±oo) = / Cj{±oo,E,£)e ^ ipj{±oo, E)e^^Q{E,e)dE. 

J A 
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With respect to (^.Sj), the only dependence left in the space variable in the integrand is in the exponent. 
The index ±00 refers to the choice of asymptotic mode A;j(iboo, E) and polarization (/7j(iboo, E) taken in the 
definition. Note the relation 

(5.10) {ied-j^y (j)j{x,t,e, ±00) = 

J Cj{±oo, E, e)e-i(fc.(±~-s)^+'^.(±oo.s))/=fcj(±oo, E)ipj{±oo, E)e-^^^''Q{E, e)dE = 

(j)f{x,t,e,±(x>). 

We also remark that since (l)j\x,t,e,±oo) are constructed as integrals in the same way as (j)^j\x,t,e) are, 
only with simpler integrands, then they also satisfy the estimates based on this structure. In particular, 
(5.6) holds without restriction on the boundary of the x-region: 

For any xq G M, there exists a constant C^{j, xq) such that for any /? G (0, 1) and any / G [0, • • • , m — 1], if 
\t\ > 1 with sign (t) = ±sign(dEkj{±oo,E)) and ±x > ±x^{j), then 



Again, assuming ( p.24| ) and ( p. 81 ), C^{j,XQ) can be chosen uniformly as e — > 0. 
Finally, (pf{x,0, e, — 00) determines (l)^j\x,0,e,+co) by means of (|3.7| ). 

While the waves (l)j{x,t,e, ±00) are not localized in space, we expect them to be approximations of 
solutions to ( [5.ID in neighborhoods oi x = ±00 only. Hence the following construction: 

Let X I— > Lo{x) £ [0, 1] be a function such that u{x) = 1 if x > 1 and io{x) = if x < 0. We define 
asymptotic waves corresponding to (j)j{x,t,e, +00) for x > 1 and to <j)j{x,t,e,— 00) for x < — 1 as follows: 

(t)f{x,t,e,a) = a;(x)(^j^'(x,t, £,+00) + (1 - w(x))0^'^(x,t,e,-cx)). 

md 

(5.12) </>W(x,t,e,a) = (x, t, a). 

i=i 

Under our hypotheses, it is easy to compute the norm of these different asymptotic states by means 
of the rescaled Fourier transform defined as: 

(5.13) {J'egmx) = [ 9{k)e-'^-l'dk. 

Lemma 5.2. Assume (HI), (H2), (H3), (CO) and (GV). Then there exists such that for all j = 
1, • • • , md, all I = 0, ■ ■ ■ ,m — 1 and all t 

\\{i£d:,y(l)j{-,t,e,±00)\\L2(JS,) < De. 

Remarks: 

i) As a direct corollary, \\{iedxy4>{-,t,e,^oo)\\ = 0(0^). Moreover, \\cl)^j\-,t,e,a)\\i2(^-^^i = 0{D^) and there- 
fore ||,/.W(-,t,e,a)||^2(u) = 0(2?,). 

ii) Again, further assuming ( |5.8D and ( p.24| ), we get = ^JeD^ with D uniform in e. 

Proof. Under (GV), the reciprocal functions oi E ^ kj{zizoo, E) all exist on A and we denote them by 
k i-^- E^{k), j = 1, • • • , md. Hence, using (|5.1[1| ) and a change of variables, we can write 



{iedx) (l)j{x,t, £,±co) = (^^■'(x,i,e,±cx)) = V27re(JP"e (•, e, ±oo))(x), 
where 

{k, t, e, ±00) = c, (±00, £:±(A:), e)e-'"^ (±~'^^*('=))/^A:Vj (±00, Ef{k))e-''''f^''^/'QiEfik),e)dkEf{k). 
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By Plancherel formula 



where is uniform in t € M. 

Finally, as expected, we show that the exact solutions (^]^) behave more and more like the correspond- 
ing free asymptotic waves ( 5.12 ), in norm, as time gets large. Furthermore, we show that ( ^.j ) cannot 
get trapped on a compact set of M as time goes to infinity, since its norm vanishes for |t| ^ oo on such 
sets: 

Proposition 5.1. Assume that (HI), (H2), (H3), (GV) and (CO) are satisfied. Then, there exists > 
such that we have for any \t\ > 0, Vj G {1, . . . , md} and V/ G {0, . . . , m — 1}; 

||4'l(-,t,e)-0W(-,t,e,a)|U2(M)<^. 



Moreover, for any bounded interval / G M, 



for some depending on I. 
Remarks: 

i) As a direct corollary. 



^"(•,i,e)llL^(/)<§, 



{ied^ycj){;t,e) - cP^^\;t,e,a)\\L2m < 



t\ 



ii) Further assuming ( ^.81 ) and ( 3.24| ), we can take = C, and Ce = C uniformly as e — > 0, see the proof. 

iii) The estimate is independent of the signs of t and of the asymptotic group velocities, because the defini- 
tion of (;(!)W(-, t, e, a) takes into account the asymptotic waves travelling in both asymptotic regions. See the 
example below for an illustration. 

iv) We prove this Proposition in 8.6. 

In order to have a better understanding of the localization properties for large times of the asymptotic 
approximation ^^'^(•, t, e, a), we need to look at the signs of the group velocities —/dEkj{±oo, E) of its 



components ( [5.91 ). Different cases occur that we list below. 

Corollary 5.1. Assume (HI), (H2), (H3), (GV) and (CO) are satisfied. Then there exists a constant Hf, 
such that and for any < (3 < 1/2 and \t\ > 1, 

{ dEkj{—oo, E)dEkj{+oo, E) < and tdEkj{+oo,E) < } 



, ,,^^\-,t,e,a) - ((/.y (•, t,e,-oo) + (•,t,e,+oo)) 11^2 
{ dEkj{-oo, E)dEkj{+oo, E) > and tdEkj{+oo,E) > } 
Uf{-,t,£,a) - (j)f{-,t,e,-oo)\\L2^^) < 
{ dEkj{—oo, E)dEkj{+oo, E) > and tdEkj{+oo,E) < } 
=^ \\(j)f{-,t,e,a) - (/>J''(-,t,e,+oo)||i2(K) < 
{ dEkj{-oo, E)dEkj{+oo, E) < and tdEkj{+oo,E) > } 
||0j^(-,t,e,a)||i2(R) < 
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t\P' 



Proof. Just make use of the definition ( |5.12| ), of ( |5.1l| ) and of the support properties of uj. 



Remark: Again, if dJ) and (1^) are true, the constant is uniform in e. 



Another consequence of Lemmas 5A and equation ( p.?]) and Proposition p.l| is the fohowing 
estimate 

Corollary 5.2. Assume (HI), (H2), (H3) and (CO) are satisfied. Then, there exists > such that for 
all I = 0, ■ ■ ■ ,m — 1, 

sup\\{edccy(J){x,t,e)\\L'2m < F^. 



If, furthermore, ( [5-4 ^ and (3.24) true, F^ can he chosen as F, uniform in e ^ 0. 

Hence, if the norm is not conserved under the time evolution (1.2), it remains uniformly bounded 
in time. Moreover, it is also uniformly bounded in e, for the type of energy densities that we will use below, 
see ( |6.1D , with { f>.8\ ) and ( |3.24D . Hence, in that case, the norm of our solutions at any time is proportional 
to that they had at any initial time to: 

(5.14) m;t,e)\\L^m) < to, e)||L2(R). 

Let us illustrate some of the notions of this Section by means of an explicitly solvable example. Consider 
the following scalar linear PDE: 

(5.15) (tanh(x)k(9t — iedx)(t){x, t, e) = 0. 

The corresponding dispersion relation yields k{x,E) = £'tanh(2;) as unique mode, which satisfies (GV). 
As the equation is e-independent, we take e = 1. The general solution reads (j){x,t, 1) = f{t + ln(cosh(x))), 
where / is any regular function. To have a solution obtained by means of a superposition of generalized 
eigenvectors e~*^^o tanh(j/)(ii/ according to some compactly supported energy density, we must have 

/(t + ln(cosh(2;))) = / e-'^^^^''''^^y)'^ye-''^^Q{E,l)dE = [ e-*^(*+''"(™"^(^)))Q(^, l)d^. 
7a J a 



Hence, / = ^/2tt{J-iQ{-, 1)) and is therefore L^, analytic and goes to zero at infinity. 

That the norm is not conserved in general under our hypotheses is now easily seen: Since < 
ln(cosh(a;)) is even and behave as |x| for x large, one checks that we have limt^_oo ||</'(-,t, 1)||l2(]{j) = 

0(||/(-)IIl2(r)) > 0, whereas, limt^+oo \\4'{-,t, 1)IIl2(r) = 0. 

Let us investigate the asymptotic waves corresponding to ( ^.151) . Using A;(ibc>o, £■) = ±E, we find 

(j){x,t,l,±oo) = [ e'^'^^e-'^^Q{E,l)dE = V2^{J^iQ{-,l)){t±x). 
J A 

Hence, 

(p{x, t, 1, a) = uj{x)V2^{TiQi-, l)){t + x) + (1 - u;{x))V2^{TiQ{-, l)){t - x), 

which, as t — > — oo, is significant at both large and positive values of x and large and negative values of 
X. Accordingly, for t ~ — oo, /(t + ln(cosh(x))) is significant at values of ln(cosh(3;)) ~ ~ i.e. for 
X ~ The picture is that of two bumps at plus and minus infinity in space that travel towards one 

another with unit velocity and disappear as they collide. This is correctly captured by the approximation 
4>{x,t, l,a), for large times. 

6. Semi-classical transitions asymptotics 

6.1. The transition integral. We assume here we are in an avoided crossing situation, and we do not 
explicit the dependence in the variable (5 > in the notation. We have obtained the asymptotics of the 
scattering matrix S{E, e) in Section ^. We now compute the small e-asymptotics of the integrals that 
describe the asymptotic states (f>j{x,t, e,±oo) given by ( |5.9D as \t\ oo, for the different channels. 



We assume that j is such that (4.1) holds and let n = '^{j) be given by (4.11). 
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We choose our energy density Q{E, e) to be more and more sharply peaked near a specific value Eq S A\(9A 
as e — > 0. As a result, we obtain semiclassical wavepackets that are well localized in phase space. This is a 
physically reasonable choice that allows for a complete semiclassical treatment. 
More precisely we consider, 

(6.1) Q{E,e) = e-^^^^l' e-'^^^'^l' P{E,e), 

where 

(CI): The real- valued function G > is in C'^(A), is independent of 5 and has a unique non-degenerate 
absolute minimum value of at Eq in the interior of A. This implies that 

G{E) = g{E-Eof/2 + 0{E - E^f , where 5 > 0. 

(C2): The real- valued function J is in C'^(A). 

(C3): The complex-valued function P{E,e) is in C^(A) and satisfies 



Qn 



dE"" 



< Cn, for n = 0, 1. 



(6.2) sup 

BgA 

e>0 

Remarks: 

i) Typical interesting choices of Q are G = g {E — Eq)"^, J = 0, and P an e-dependent multiple (the equation 
(1.2) is linear) of a smooth function with at most polynomial growth in {E — Eq)/£. 

ii) We want to emphasize the fact that a Gaussian energy density does not give rise in general to a Gaussian 
solution. See the discussion in the Introduction and Section 6 of IH 



The leading inter-modes transitions are described by the asymptotics of those coefficients q(iboo, E^ e) that 
satisfy 

(6.3) Cfc(-oo,^,e) = 5j^k 

(6.4) Cn{+oo,E,e) = e'^^^^^'^) e'-^c (1 + 0£;(e)), 

where n = 7r(j) = j it 1. We recall that the error term Oe{£) depends analytically on the energy £^ in a 
neighborhood of the compact set A. We have already noted in the comments after Theorem iA that the 
term Oe{£) satisfies (|6.2[ ). 



Theorem 6.1. Assume (HI) to (H3), (AC) and (GV). LetQ{-,e) be the energy density supported on the 
interval A defined in l\6. 1\) which satisfies (CI), (C2), and (C3). Let (j){x,t,£) be a solution of equation 



(5.1) of the form (5.i). Assume dEkj{— 00, E) < on A, for some j and suppose that the solution is 
characterized in the past by 

lim \\(t){-,t,e) -(/){-, t,e,a)\\L2(^) = 0, 

where, as t ^ —00, 

ij{x)(j){x,t,e,a) = [ Q{E,e)e-'*'^/'e-'^'^''^''-°^'^^+'^^^-°°'^^^/^ipj{-oo,E)dE + Oil/\tf). 
Let n = TT(j) be given by ( 4-11 ), and let 

(6.5) a{E) = G{E) + lm{j^kj{z,E)dz), 

(6.6) k{E) = J{E) - Re{j^kjiz,E)dz) + uJn{+oo,E). 

Assume E* is the unique absolute minimum of a{-) in Int A and define k* = kn{+oo, E*). Let k 1-^ E^{k) 
be the inverse function of E ^ kn{±oo,E) on A. 

Then, there exist 6q > 0, p > arbitrarily close to 5/4, and a function eq : (0,(5o) such that for all 

< (3 < 1/2, 6 < 6q, and e < Sq{5), the following asymptotics hold as t ^ —sigw. {d^E^ {k*))<x , in L"^' 
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norm: 



(6.7) <pn{x,t,e) = ^/2;^P(^*,e)e-"(^*)/^e-^'^(^*)/Vn(+oo,^*)e-^^^(f'^*)9fc£;+(r 



where 



A{k) = y (A: - k*f + iXi{k - k*), with Ai = dkE+{k*)K'{E*), 



A2 = [4i?+(fc*)]'a"(i?*) +i [^"(i?*)[4i?+(A:*)]2 + K'(i^*)a2£;+(r)] , 

c^^fi^ Xk„{+oo,A) is the characteristic function of the set kn{+oo, A). 
Moreover, if t ^ sign {dkE+{k* ))oo, then \\(pnix,t,e)\\L2^^) =0(l/|t|^). 

Remarks: 

0) The first error term is uniform in t whereas the second error term is uniform in e. 

1) The same result holds for (f>, and 4>n{+oo), replaced by (iedYcj), and {iedY (f>n{+oo) respectively, with 
{I = 0, - ■ ■ ,m — 1}, at the expense of a multiplication of the prefactor by kn{+oo, E*y . 

a) As will be made explicit in Section |^ below, the norm of the leading term expressed as a Fourier 
transform is positive, of order e^^^, and independent of time. The leading term hence becomes meaningful 
for times t that are of order |t| ~ e"^/^, for some c > at least. We get control over this time scale far beyond 
the Ehrenfest or Heisenberg times of Quantum semiclassical analysis thanks to our scattering setup. 

iii) The leading term clearly satisfies the asymptotic PDE ( 1.14 ). 

iv) The energy E* depends explicitly on the properties of the involved modes and on the energy density 
Q{E, e) as well. 

v) The space-time localization properties of the leading term are further discussed in Section |^ 

vi) Also, as mentioned earlier, we can specify the coefficients q at x = +00 instead. 

vii) The proof of the Theorem is given in the last Section of the paper. 



Let us finally discuss our hypotheses and interpret our result. The condition on the sign of dE{kj (—00, E)) 
says that the group velocity of (f)j{x,t^e,— 00) is positive, so that <j){x,t,£,a) is non trivial as t ^ —00 for 
negative x's and describes an ingoing wave. If the asymptotic group velocity of the mode kn{x, E) is positive 
as X — > +00, our results describes an outgoing transmitted wave for large positive times, as discussed in the 
introduction. If the asymptotic group velocity —dkE^{k) is negative, we describe another ingoing wave 
along mode n, for large negative times and large positive x's, arising during the evolution, which, as time 
goes to +00, goes to zero. Note also that if the asymptotic group velocity of mode kj{x,E) at a; = +00 is 
positive, then an order one wave, in the sense that Cj{+oo,E,e) = 1 + 0(e), propagates along positive x's 
for positive times. If the asymptotic group velocity of mode kj{x,E) at x = +00 is negative, there is no 
wave propagating along positive x's to the right, for large positive times, but another ingoing wave from 
large positive x's and large negative times. 

Therefore, in case —dkE^{k) is negative and both —dkEj{k) and —dkEj~{k) are positive, running 
the evolution backwards in time, we have an ingoing wave (of order one in the sense above) on mode j, 
for X — > +00 and t — > +00, and, as t — > —00, we have an outgoing wave on mode j, for x — > —00, and 
another exponentially small outgoing wave on mode n, for x +00, whose asymptotics is determined by 
our Theorem. Hence, we describe the asymptotics of a reflected wave in mode n. Note that reflected waves 
on other modes may be present as well. In any case, they are exponentially small. 

Finally, in case —dkE^{k) and —dkEj~{k) are both negative, we describe a scattering process in which 
we have ingoing solutions on the modes j and n, that all disappear as time goes to +00, in a similar way as 
what happens in the illustration ending the previous Section. 

6.2. Perturbative results in 6. We assume that (H4) is also satisfied and restore back 6 in the notation. 
We have the following sharper result concerning the behavior as 5 —> of the quantities involved in the 
description of the asymptotic wave: 
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Proposition 6.1. Further assuming (H4), we have the following as 5 ^ 0, for E G A; 



Im j{ki - kj){z, E, 6) dz = D{E)6^ + 0{6^), 

wuh D{E) = r'^^^'liflf^^\ 

This implies that (E, 5) ^ Im J^{ki — kj)(z, E, 6)dz is a positive function. 
Let a{E,6) = G{E) + Im J^{ki — kj){E,5). There exists E*{6) such that 

dEa{E*{6),6) = 0. 

It satisfies: 

E*(5) = Eo-^^^6^ + Oi5^). 
9 

The results above hold provided one knows E*{5) is the unique absolute minimum of a in the set A, 
which is generically true. Again, if there are several minima, one simply adds the corresponding contri- 
butions. Note also that if the constant g characteristic of the energy density is of order 5^, the difference 
Eq — E* is of order one as 6 shrinks to zero. This corresponds to a "wide" energy density of width 
around Eq. This result is a straightforward consequence of the Implicit Function Theorem, the proof of 
which we omit. 



6.3. Explicit computation in case E^{k) is quadratic. In this paragraph, we assume that k i— > E^[k) 
is quadratic: 

(6.8) VA: G A;„(A, +cx)), dlE+{k) = 0. 

This is true for all modes in the study of the Born-Oppenheimer approximation, see |19]. This situation 
allows for an explicit determination of the leading term in the asymptotic wave. We also assume that the 
function a has a unique absolute minimum E*{5). For sufficiently small 5, this minimum is non degenerate 
and satisfies E*{5) G IntA. 

The following result is proven in Section 

Lemma 6.1. Assume that k E^{k) is quadratic and that a has a unique absolute minimum E*{5) G /niA. 
There exists p G]3/4, 5/4[ such that, as e ^ Q: 



(6.9) 0„(x,t,e,+oo) = e-"(-^*)/^e-^'^(^*)/V„(+cxD,£;*)e-*^^(f'^*)p(^*,e)afcK+(r) x 



[\2 + idlEt{k*)tYl'' 



e 



Remarks: 

i) The leading term in that case is a freely propagating Gaussian, i.e. an exact solution to 



iedtg{x,t,e) = (^E* + dkE^{k*){ied, - k*) + 



^^^ized,-kr]9ix,t,e), 



centered at Xc{t) = —dkE^{k*)t — Ai, of width y/ei and of Lp' norm of order e^/^. 

ii) In the general case, the error terms involved in the course of the computation are not uniform in time, 
which prevents us to get such an explicit form for the asymptotic wave. Nevertheless, we show in the next 
Section that we can get a fairly accurate description of such asymptotic waves, for large times and small e. 
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7. Space-time properties of the asymptotic waves 



As seen above, the interpretation of our results makes use of the space-time properties the different 
asymptotic waves (l)j{x,t,e, ±00) in terms of which the time-dependent scattering processes are expressed. 
The present Section is devoted to a thorough description of the space-time properties of the leading term of 
these waves as e ^ and Itl — > 00. 



We first note that Theorem |6.1| also holds for the wave t, e, ±00), which characterized by the 
asymptotics Cj{±oo, E,e) = 1 -|- Oe{£)- It suffices to replace the index n by j, the values E* and k* by Eq 
and ko, and to set a{E) = G{E) and ^j(Ci E) = 0. Note in particular, that a(£'o) = 0, as it should be. 

Therefore, the space-time properties of the asymptotic waves along modes j and 7r(j) = n are encoded 
in the Fourier transform 

1 



(7.1) 



-i{kx+tEf{k))/e^-K{k)/e 



27re Jki{(7,A) 

where the index I stands for j or n, and a for -|- or — . We will also denote k* or k^, respectively E* or 
£"0, depending on the context, by k, respectively E. We can make use of the positivity of the real part 
of the function A(A;) and of Parseval's formula to regularize and localize the integrand as follows. Let 
rj G Cq°(]R) with support in [—1, 1] and r]{k) = 1 in a neighborhood of /c = 0. Set r]i;{k) := r]{{k — k)/e'^), 
with < r < 1/2. Then, if k ^ supp(%), le-^C')/^] = 0(e~). Therefore, we have in L norm, 

1 

Xki (+00,. 



-iikx+tE^{k))/e^-A{k)/e^^^j^^ dfc + 0(£ 



(7.2) 



where the error term is uniform in t. Note also that by Parseval again 
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(7.3) 




dz + 0(e= 



[dkEr{k)]^a^^{E) 



+ 0(e°°). 



uniform in t. Hence, the norm of the asymptotic state in Theorem |6.1| is positive, independent of time 
and of order e^/^ . 

Now, as ki{aoo, •) is analytic m E € A, the same is true for the inverse function -Bf (•) in A; G ki{+oo, A). 
Moreover, e~^^'^^^r]^(-) is in Cq°, so that we can apply stationary phase methods to describe the large t and 
X behavior of ([7^. 



Proposition 7.1. Let rji; be as above and 1 > a > 1/2 and assume dkEf{k) 7^ 0. Define for all \t\ > 1, 



Ct{e) 



U 



k~k\ <e 



^{x £R\\x + dkE[{k)t\ < 



Then, there exist > and c{n) > 0, such that for all e < Eq, all n and all \t\ > 



(7.4) 



/27re 



-i{kx+tEr{k))/e^-A(k)/e^^^j^^^j^ 



L2(R\Ct{e)) 



< c(n)- 



O 



t|3a/2-l 

e 



e 



\2a-l 



2a- 1 



Remarks: 

i) The Proposition says, essentially, that the whole L?' mass of the asymptotic wave in Theorem |6.1| is located 
at time t in (a slightly larger) neighborhood of size \/\t\ of the point propagating with the group velocity 
—dkE^{k), up to arbitrarily small corrections as £/\t\'^'^~^ 0. 

ii) The Proposition actually also holds if e = 1, if one is not interested in the small e behavior. 

iii) The condition |t| > l/e^^^^~"^ actually represents no restriction in our case, since we need to work with 
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exponentially large times in e, in order to have a meaningful leading order term in Theorem 6J. 
iv) The proof is given in the last Section. 

While we don't need to assume anything on the direction of propagation of the involved waves for 



Theorem 3.1 to hold, it's usefulness in describing time-dependent scattering processes is revealed by the 



above interpretation based on these directions of propagation. 

8. Technicalities 



8.1. Proof of Lemma |2.1j . We first prove that if A is small enough, p is independent of E. Fix Eq ^ /\ 

and xq such that {ki — kj){xQ, Eq) = 0. By hypothesis ki — kj is continuous, then, by Cauchy formula, 
dxiki — kj) is continuous and dx{ki — kj) ^ in a neighborhood of {xo,Eq). By local inversion, the set 
{E G A ; p{E) = p{Eq)} is open. Thus, for any E £ A, p{E) is constant. 

By linear perturbation theory, there exists y > such that there are no non real crossings in /jy, for any 
Ee A. 

We have the following result: 

Lemma 8.1. Fix i ^ j. The functions ki and kj have the following properties 

The function {z,E) t-^ ki{z,E,0) + kj{z,E,0) is analytic on py x A. 
(2) The function {z,E) ^ [ki{z, E,0) — kj(z,E,0))^ is analytic on py x A. 

Proof. According to [^6|, we know that we only have to check the analyticity of ki + kj and {ki — kj)"^ in 
a neighborhood of a crossing, actually a branch point {xq,Eq). Let P{z,E) be the 2-dimensional projector 
on the A-group corresponding to the eigenvalues ki{z,E,0) and kj{z,E,0). Let T be a small a closed path 
in C surrounding ki{xo, Eq). For (z, E) in a neighborhood of (xq, Eq), we can write : 



P{z, E) = ^ l^{H{z, E, 0) - X)-'dX. 



As H and thus its resolvent on F are analytic in (z,E), this implies that {z,E) i-^ P(z,E) is analytic in a 
neighborhood of {xq, Eq). We consider {ipi{xQ, Eq), (p2{xQ, Eq)} a basis of P{xq, Eq)C"^'^ and we define: 

ipi{z,E) = P{z,E)ipi{xQ,EQ) ; ip2{z, E) = P{z, E)lp2{xq, Eq). 

Then, in a neighborhood of {xq,Eq), {ipi{z, E), ip2{z, E)} is an analytic basis of P{z, E)C"^'^. The matrix 
M{z,E) oi P{z, E)H{z, E,0)^p(^^^^-^Qmd on the basis {ifi{z,E),ip2{z,E)} has analytic coefficients. Besides, 
a{P{z, E)H{z, E, 0) |p(^^£;)C'»d ) = {h{z, E, 0) , kj {z, E,0)}. This implies that det M{z, E) = ki{z, E, 0)kj {z, E, 0) 
andti M{z,E) = ki{z, E,0) + kj{z, E,0). We finish the proof with the identity (ki-kj)'^ = {ki + kj)"^ — Akikj. 
This ends the proof of Lemma |8.1| . 

We define: 

A = {{z,E)/3i / j7 {ki - kj){z,E,0) = 0} C M x A. 

According to |26| |, it suffices to prove that the functions {kj}j=i^,,,^rnd are analytic in a neighborhood of 
any (xo,-E'o) G A. Fix {xq,Eq) £ A. There exist i and j such that ki{xQ, Eq,0) = kj{xQ, Eq,0). For I / j 
and I ^ i, ki is analytic in a neighborhood of {xq,Eq). By using Lemma |8.l| , it suffices to prove that 
{z,E) {ki{z,E,0) — kj{z, E,0)) is analytic at {xq,Eq). The function g{z,E) = {ki — kj)'^{z,E,0) is 
analytic in a neighborhood V of {xq, Eq). Besides, since ki — kj is real for any {x, E) in y n M^, we have: 

g{x,E)>0, V(x,^)GynM^ 

We write the Taylor expansion of {z, E) g{z, E). There exist (a, f3, 7) G such that: 

g{z, E) = [a{z - xq)^ + - xq) + j{E - Eq)]^ + o{{z - xq)^ + {E - Eq)^). 

We start with proving that there exist 7) G M? such that 

g{z, E) = [p{z - xq) + ^{E - Eq)]^ + o{{z - xq)^ + {E - Eq)^). 
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We know that, for any E £ A, the function z ^ \/giz, E) is analytic. For \E — Eq\ and \z — xq\ small 
enough, we have that: 

The function z t-^ \/ g{z, E) can be analytically continued in a neighborhood of xq only if 7a = 0. This 
proves the announced result, with 7) = (/?, 7) or (/?, 7) = (a/q;^ + Z^^, 0). 

We notice that we have the following relations: 

\d^{ki-k,){x^,Eo)\ = \P\^Q. 

\dE{ki - kj){xQ,EQ)\ = |7|. 
To end the proof, it remains to show that: 

g{z, E) = [P{z - xo) + j{E - Eo)]^ + 0(|/3(z - xq) + j{E - Eo)f + \E- E^Mz - xq) + 7{E - Eo)f). 

We change variables for u = j3{z — xq) + '^{E — Eq) and e = {E — Eq). 
Since P ^ 0, this map is bijective and we consider the function g: 

~i \ fu-^e ^ \ 

g[u,e) = g y — ^ '^xq^c + Eq j . 

We write the Taylor expansion of g near (0,0): 

g{u,e) = + ^ aiqU^e''. 

l+q>3 

Since g is real positive on a neighborhood of (0,0), we obtain that aiq = 0, for any q>2. This implies: 

g{u,e)=u^ll+ aiqU^^'^eA =u^{l + 0{\e\ + \ue\)). 

\ '>2,(J>1 / 



Point (3) is an immediate consequence of (2) and of Rouche's Theorem. It ends the proof of Lemma 2.1 



8.2. Proof of Lemma |2.2| . A proof similar to the proof of Lemma |8.1| shows that the function (z, E) 
iki — kj)'^(z, E,6) is analytic for 5 small enough and that {z,E,6) 1— > {ki — kj)'^(z, E,6) is C^. 



For 5 = 0, Lemma 2T implies that a{E) = \dz{ki — kj){x(){E), E,0)\ is analytic in E £ A. For 5 > 0, we 
define: 

ih - kjfiz,E,5) - a\E){z - xo{E))^ 



ri{z,E,6) 



25 



The function E 1— > dzri{xQ(E), E, 6) = c{E) + 0{5) is analytic and 0{6) is uniform in E. Thus we can apply 
Weierstrass's Theorem to get that the function c satisfies c{E) = lim^^o dzri{xo{E), E, 6) and is analytic in 
E e A. We also define: 

ih - kj)\z,E,S) - a\E){z - xo{E)f - 2c{E){z - xo{E))6 
r2{z,E,d) = ^2 • 

Similarly, the function 6^ such that h'^{E) = Wm^^Q r2{xQ{E) , E , 5) is analytic. 
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8.3. Proof of Lemma |3.3| . We fix 5 > and drop it from the notation. For p > 0, we denote by 
Vp{A) = {£ £C; dist(£:, A) < p}. 

By perturbation theory, there exists Y > and p > 0, depending on 6, such that {z,E) ^ kj{z,E) is 
analytic on py x Vp{A) and: 

inf ihiz.E) - kJz,E)) > 0. 

t^j,{z,E)&pYxV,{Ay 

Hence, for a > small enough, we can write the eigenprojector Pj{z,E) as: 

Pj{z,E) = ^ [ [Hiz,E)-X]~'dX. 

J\X-kj{z,E)\=a 

We recall the identity 

[H{z,E) - A]~^ - [H{oo,E) - X]-'^ = [H{z,E) - X]~^ [H {oo, E) - H{z, E)][H {oo, E) - A]"^ 

and the fact that H{-,E) and its derivatives with respect to E satisfy the analog of H2 uniformly in 
E E V^(A). This implies that for any j G {1, . . . , md}, for any I £ N, and uniformly in E £ Vp{A): 

(8.1) sup |Rez|2+^||a^(Pj(z,S) -Pj(+c30,S))|| + sup |Rezp+''||aij(Pj(z, ^) - Pj(-cx), ^))|| < cx). 

Rez>0 Rc2<0 

We consider the vector $j(0, E) satisfying ( ^.141) for z = and which is analytic in E. Consider the identities 
k,{z,E)=tT{Pj{z,E)H{z,E)), <l>j{z,E)=W{z,E)^j{0,E), 

and Cauchy formula: 

for r and r' small enough. Then, ( 3.171) fo llows. To get ( 3.1S| ), one also makes use of the differential equation 
satisfied by W{-,E), of the estimate ( |8.1j ) and one considers the first d components of ^j{z,E) only. 

The explicit form ula (|3.9D and the decay of dxkj and of dx^j yields formula (|3.2[l| ). 

Finally, estimate ( |3.23| ) is a direct consequence of the definitions (|3^) and ( ^.22^ and of ( p7|) . 

8.4. Proof of Lemma |3.4| . The proof of Lemma 3.4 is virtually identical to the one of Lemma 3.1 in |19|, 
once the properties of the matrix M{x,E,£) have been established. Therefore, we just give here the main 
steps of the first part of the argument for the reader's convenience. 

Because of ( |3.2Cl| ), we know that there exists C such that, uniformly in E and e, 

\\M{y,E,e)\\dy <C. 
Expressing the solutions of equation ( |3.7| ) as Dyson series: 

.2) c{x,E,e) = 



JO 



n=0 ■ 

X M{xi,E, e)M{x2,E, e) ■ ■ ■ M{xn, E, e)dxidx2 ■ ■ ■ dxn c(0, E, e) 

we obtain the usual bound: 

\\c{x,E,e)\\ < ei'<ni*^(^'^'^)ll'^^ ||c(0,£;,e)||. 



Thus, we get from ( ^.2D that c{x,E,£) is bounded as x — > ±oo. Next we show that \\c{x,E,£) — c{y,E,£)\\ 
is arbitrarily small for large x and y, so that limx^cx)C{x,E,e) = c(oo,£', e) exists. It is enough to consider 

ry 

c{x,E,e) — c{y,E,e) = — M{z, E,e) c{z, E,e) dz. 

J X 
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The expression above with y = ±00, and the properties of M, c, just proven yield the bound 



c{x, E, e) - c(±oo, E, e) = 0{< x >-^^+''^). 
Finally, if ||c(iboo, e)|| is further assumed to be uniformly bounded in i? G A and e 



0, it is enough to 



consider the initial conditions (4^), by linearity. Then, by integration by parts, see (|4.7| ) and (|4.8| ), we get, 

c{x,E,e) = 0(1), 

uniformly hi E £ A, e ^ and x G M. Hence, all bounds above are uniform in £' G A and e ^ 0. 

Getting similar bounds on the derivatives of c with respect to E which are uniform in e and E £ A requires 
a little more work. The argument is identical to that used in |19]. We resort again to integration by parts 
in ([4. 61) with xq = —00, differentiate with respect to E, and make use of Gronwall Lemma to get bounds. 



We do not give the details and refer to Section 7 of |19]. 



8.5. Proof of Lemma 5.1. Again, we mimic the proof of Lemma 4.1 in It suffices to rewrite: 



dg-i{Jo kjiy,E)dy+tE)/e 



l€ 



We compute: 

cl)j{x,t,e) 



IE 



le 



+ ie 



Cjjx, E, e)Q{E, e)^j jx, E, e) kj{y,E)dy+tE)/, 
t + Jq dEkj{y,E)dy 

dEjCjjx, E, e)Q{E, e)ipj {x, E, e)] ^_^( k,{y,E)dy+tE)/e^^ 

t + dEkj{y,E)dy 
Cj{x, E, e)Q{E, e)ipj (x, E, e)e-'^Io kj{y,E)dy+tE)/. 



[t + f^dEkj{y,E)dyY^ 



df^kj{y,E)dydE. 



By Lemma 3.3 we have for x large enough: 

r r t /I 

t+ dEkj{y,E)dy = x dEkj{±oo, E) + - + O { - 
Jo I X \x 

Under the restrictions put on x/t, we obtain that for a constant independent of t: 



ii^,t,e)\\ < 



X 



Now, for any / G {0, • • • , m — 1}, we have: 

rn f , -i(J^ kj(y,E)dy + Et) 

cl)y{x,t,e) = J Cj{x,E,e)k){x,E)Lpj{x,E)e ^ Q{E,e) dE, 

so that a similar computation gives the result. 

The last estimate makes use of the fact that under the given conditions stated on the signs of x and t, 
and for \x\ large enough, uniformly in time and energy. 



(8.4) 



\t+ [ dEkjiy, E)dy\ >\t + xdEkj{±oo, E)\/2 > {\t\ + \x\)c, 
Jo 



for some c > 0. The result follows from the elementary inequality (x + y) > x^y^ ^, for any x,y > and 
any f3 G (0, 1), and the arguments used above. 
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8.6. Proof of Proposition 5^. We adapt the proof of Proposition 4.1 in |19| and we only give the main 
steps. 

i) We rewrite: 



.5) </.W(x,t,e)-0^'J(x,t,e,±oo) 



[I] I 



[k]{x,E)ipj{x,E) - k\{±oo,E)ipj{±oo,E)]e 



'i So kj(y,E)dy 



kU±oo, E)ipj {±oo, E)e 



-i(kj(±co,E)x+LO^(±ca,E)) ■ ±r t^-, i 

— 7 — fe~*''j ^^'^11^ 



Cj{x, E, e)Q{E, e)e-**^/^ dE + 
l]cj {x, E, e)Q{E, e)e-'^^/^ dE + 



A;j(±oo, E)^pj{±oo, E)e' 



■ [cj{x, E, e) - Cj (±cx), E, e)]Q{E, e)e-^*^/" dE. 



ii) To obtain a factor 1/t, we integrate by parts. For any regular function /, we have: 



(8.6) 



f{x,E,e)e-'^^/^dE 



l£ 



f{x,E,e) e-**^/^ 



dA 



le 

T 



dEf{x,E,e) e-'^^/'dE. 



iii) We apply equation ( ^.61) for the three terms in 



iv) The estimates of Lemma 3.3 and Lemma |3.4| prove that there exists Cg > such that: 

a 



sup \x\''\\[4>f{x,t,e) - 4>f{x,t,e,±^)]\\ <^ 

X— >±oo ll' 



v) The identity {x,t,£) = uj{x)(l)j {x,t,£) + (1 — uj{x))(j)j {x,t,e) and support considerations on the 
definition of 4>^j\x,t,e,a) yield the first result. 

vi) The second estimate follows from ( |8.6D applied to the definition of (j)^j\x,t,£). 

8.7. Proof of Theorem OA. Taking into account the results of Section ^, we are left with the computation 
of the following integral: 

J A 

with P{E,£) = P{E,£)e~^^^^'^'^\l + Oe{£)), where 0£;(e) is defined in (6^). We already remarked that 
P{E,e) satisfies C3 as well. Again, wa adapt the arguments of [p^]. 

In terms of the variable k = kn{+oo, E), we rewrite T as: 

ife„{A,+oo) 

We can see T as the rescaled Fourier transform (see ( |5.13| )) of the function: 

S{£,k,t) = v/2^P(^+(A;),e)e-"(^"(+~''=))/"e-^^(^"('=»/^e-^*^"('=)/Vn(+oo,^+ 
We follow the analysis done in ||l9| and expand S around k* 



a{Etm = a{E*) + [dkE^ {k*)fa^E*) ^^ + 0{{k - k*f), 



where a"{E*) > 0. We define: 



K{k) = K{E+{k)), 



k{k) = k{k*) + k'{k*){k - k*) + K^k*) ^^ + 0{{k - k*f). 

We also have: 

P{E+{k),E) = P(i^+(r),e)e-*^^(<'^" + 0{{k - k*)) + 0(e). 

(afc£;+)(A;K(+oo,£;+(A;)) = (afcE+)(r)<^„(+oo,ii;+(r)) + o((fc-r)). 

We can then follow the computations of H^. We set //(e) = e*, with l/3<s<l/2, so that ^{ef" /e « 1 << 
f^{e)/^/£, as e — > 0. Using Parseval formula and thanks to the properties of a, the restriction of the integration 
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set to the interval [k* — k* + /x(e)] causes an error, in the sense, of order e~'^i^*)/'^£°° , On that set, 
the overall error err(e) stemming from the expansion around k* is err(e) = 0{e + + = 0(e^*~^), 

with our values of s. The norm of the error term can thus be estimated by Parseval formula to yield 



erri 



uniformly in t. The interval of integration can then be restored to its initial value kn{+oo, A) at the expense 
of another error of order e~"(^*)/^e°°. Therefore, we obtain that there exists 3/4 < p < 5/4 so that in 
L^-norm and with an error term that is uniform in t: 



(8.7) 



T{e,x,t) = e 

X 



-A{k)/e^-itE+{k)/e^-ikx/e^j^ _|_ Q^^~a(E'')/e^py 



fc„(+oo,A) 



8.8. Proof of Lemma 6.1. To prove Lemma 3.1, it suffices to compute: 

^-Aik)/e^-itE+(k)/e^-ikx/e^j^_ 

fcn(+oo,A) 



According to (6.8), we have: 



Et{k) = E+ik*) + d,E+ik*)ik - k*) + M^^(A; - k*f. 



With the same arguments as in [|19[, we can expand the interval of integration to the whole line at the 
expense of an error that is uniform in t, and in the sense: 



/ 



^~A{k)/e^-itE+{k)/e^-ikx/e^j^ 

I k„{+co,A) 

The rest of the proof follows from the formula: 



+ 00 



-A{k)/e^-itEi{k)/e^-ikx/e^j^ _|_ 0{e°°). 



with 



-(M(k-k*f /2+iN(k-k*))/e ^ 

M = \2 + idlE+{k*)t, 
N = Xi + dkE+{k*)t + x. 



e2TT 



8.9. Proof of Proposition 7.1. We will make use of the following Lemma whose proof we omit. 
Lemma 8.2. Let f and be C^(M, M). Provided f'{k) ^ on supp Us, we have for any n G N 

'u,{k)\' 1 V 



~'^^'''^/'u,{k) dk = {-ieT / e-'^^^^/' 



mj f'{k) 



1 



f'ik) 



dk , 



where there are n factors 1/ f'[k). Moreover, setting g = 1//', there exist constants c{j,ri, ■ ■ ■ ,r„) G R such 
that 



{{Ue9)' 9 



9 



1=0 



where h^^"^ denotes the k^^ derivative of h 



In our case, f'{k) = x + dkE^ {k)t, so that explicit computations show the existence of smooth functions 
on Ss, the support of = e~^('^)/'=r/£, k Cj^n{k) which are independent of x and t and such that 

(n) n 
j=0 
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1 



{x + dkEf{k)ty+^- 



Then, we get from the above 



Lemma 8.3. Assume dkE^{k) ^ and consider {x,t) such that inf^g^^ \x + d^E^ {k)t\ > 0. For any n G N, 
there exists a positive constant C{n), uniform in {x,t) and e such that 



(8.9) 



„ n 
Jm. n 



my 



(inffces^ \x + dkE[{k)t\) 



r+n ' 



for all £ < £q, where £q is uniform in n and {x,t). 



Remark: The Lemma actuahy holds for any value of e if one is not interested in the behavior in that 
parameter. 

Proof. One first notes that by scaling and by our choice of r the successive derivatives of 

satisfy sup^gj^ u\^\k) = 0(l/e-^), if e is small enough. Then, formula (|8.8|) and the restriction Y^^=i "^i — 3 ™ 



the second point of Lemma |8^ give the possible numbers of factors (^x+dkE" {k)t)i+'^ ™ integrand. Finally, 
the fact that all functions of k involved are smooth, together with \Sf;\ = e'^ yield the result. 

To end the proof of the Proposition, we note that on the complement of Ct{e), the right hand side of 
( |8.9D is bounded above by 

C(n)e^e|t|i-° / e ^ ""^ 



imikesAx + dkE[ik)t\) \\t\^^-\ 

using the fact that elil""^"" > 1. Explicit computation of the norm of (inffcg5^ \x + dkE[ {k)t\) on the 
complement of C((e) yields the result. 
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